Chiral two-body currents in nuclei: 
Gamow- Teller transitions and neutrinoless double-beta decay 
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We show that chiral effective field theory (EFT) two-body currents provide important contribu- 
tions to the quenching of low-momentum-transfer Gamow- Teller transitions, and use chiral EFT 
to predict the momentum-transfer dependence that is probed in neutrino-less double-beta (QvfiP) 
decay. We then calculate for the first time the Of/3/3 decay operator based on chiral EFT currents 
and study the nuclear matrix elements at successive orders. The contributions from chiral two-body 
currents are significant and should be included in all calculations. 
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Weak interaction processes provide unique probes of 
the physics of nuclei and fundamental symmetries, and 
play a central role in astrophysics The structure of 
strongly interacting systems is explored with /3 decays 
and weak transitions. Superallowed decays allow high 
precision tests of the Standard Model, and 0^/3/3 decays 
probe the nature of neutrinos, their hierarchy and mass. 
Weak processes mediate nuclear reactions that drive stel- 
lar evolution, supernovae and nucleosynthesis. 

Surprisingly, key aspects of well-known decays remain 
a puzzle. In particular, when calculations of Gamow- 
Teller (GT) transitions of the spin-isospin-lowering op- 
erator gAcrr~ are confronted with experiment, some de- 
gree of renormalization, or "quenching" q, of the axial 
coupling = qgA is needed. Compared to the single- 
nucleon value qa = 1.2695(29), the GT term seems to be 
weaker in nuclei. This was first conjectured in studies of 
/3-decay rates, with a typical q sa 0.75 in shell-model (SM) 
calculations [2J and other many-body approaches [3j- In 
view of the significant effect on weak reaction rates, it is 
no surprise that this suppression has been the target of 
many theoretical works. It is also a major uncertainty 
for 0^/3/3 decay nuclear matrix elements (NMEs), which 
probe larger momentum transfers of order the pion mass, 
p ~ rrin, where the renormalization could be different. 
Here we revisit this puzzle based on chiral EFT currents. 

Chiral EFT provides a systematic basis for nuclear 
forces and consistent electroweak currents 0, Hj], where 
pion couplings contribute both to the electroweak ax- 
ial current and to nuclear interactions. This is already 
seen at leading order: qa determines the axial one-body 
(lb) current and the one-pion-exchange potential. Two- 
body (2b) currents, also known as meson-exchange cur- 
rents, enter at higher order, just like three-nucleon (3N) 
forces As shown in Fig.[TJ the leading axial contribu- 
tions are due to long-range one-pion-exchange and short- 
range parts with couplings 03,04 and cd, which also 
enter the leading 3N (and subleading NN) forces [1, Q . 
Although the importance of 2b currents is known from 
phenomenological studies , chiral currents and the con- 
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FIG. 1: Chiral 2b currents and 3N force contributions. 

sistency with nuclear forces have only been explored in 
the lightest nuclei [H, H, ■ In this Letter, we present first 
calculations for GT transitions and for the CV/3/3 decay 
operator based on chiral EFT currents. A preview of the 
NMEs (Fig. [2|) and the quenching of gA (Fig. [3|) shows 
the great importance of using chiral 2b currents in nuclei. 

At low energies, the coupling to weak probes 
is given by the current-current interaction, Hyy — 
^| / d 3 r eT %v ' v ihiiJ^} +h.c, where Gf is the Fermi con- 
stant, p the momentum transferred from nucleons to lep- 
tons, and the leptonic current of an electron coupled 
to a left-handed electron neutrino [l|. In chiral EFT, the 
nuclear current J^J is organized in an expansion in pow- 
ers of momentum Q ~ m^ over a breakdown scale Ab ~ 
500 MeV. Consistently with nuclear forces [4J , we count 
the nucleon mass as a large scale, corresponding numeri- 
cally to Q/m ~ (Q/Ab) 2 , so that the leading relativistic 
1/m corrections are of order Q 2 , and 1/m 2 terms of order 
Q A . To order Q 2 (and also Q 3 in this counting), the lb 



current, jf (r) = £f=i V [6»° J° 
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has temporal and spatial parts in momentum space [5j 
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where E — Ei — E[, p = p, — p^, and P = p^ + p^; and 
vector (V), axial (A), pseudo-scalar (P), and magnetic 
(M) couplings, gv(p 2 ), 9a(p 2 ), 9p(p 2 ), and g M (p 2 ) M- 
In chiral EFT, the p dependence is due to loop corrections 
and pion propagators, to order Q 2 : 9v,a(p 2 ) = 9v,a (1 — 
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FIG. 2: (Color online) Nuclear matrix elements M 0vPI3 for 
Ov/3/3 decay. At order Q°, the NMEs include only the lead- 
ing p = axial and vector lb currents. At the next order, 
all Q 2 lb-current contributions not suppressed by parity are 
taken into account. At order Q s , the thick bars are pre- 
dicted from the long-range parts of 2b currents (cz> = 0). 
The thin bars estimate the theoretical uncertainty from the 
short-range coupling cd by taking an extreme range for the 
quenching (see text). For comparison, we show the SM re- 
sults of Ref. [1 21 ] based on phenomenological lb currents. 
The inset (representative for 136 Xe) shows that the GT part, 
M gt " =/*C'gt(p), is dominated by p ~ 100 MeV. 
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m n = 138.04 MeV, and g^ pn = 13.05 [Tl|. 

At leading order Q , only the momentum-independent 
gA and gv terms contribute. They give rise to p < 
lMeV GT and Fermi (t~) single- /3 and 2vf3(3 decay. 
On the other hand, when studying processes that probe 
larger momentum transfers, such as Of/3/3 decay with 
p ~ 100 MeV, terms of order Q 2 need to be included. For 
01/(9/3 decay, the Q° terms are still most important and 
the axial term dominates. In SM calculations, one has 



M, 



>° .axial 



1.20 and M\ 



0.15 



>P , vector/ 

3|. Other calculations 



compared to the final M 0vt313 [T 
yield similar ratios, also for the Q 2 terms, except the 
vector term can be significantly larger [l3l . [l4| . 

Among the Q 2 terms, form-factor-type (ff) contribu- 
tions and the gp part of J^ib dominate. The pseu- 
doscalar term is important, because pgp{p 2 ) ~ 7.9 for 
p ~ 100 MeV in 0^/3/3 decay. They reduce the NMEs: 
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The remaining Q 2 terms are odd under parity, so they re- 
quire either a P-wave electron (whose phase space is sup- 
pressed by 0.03 - 0.06 for Oi//3/3 decay candidates) 
or another odd-parity term to connect + states. There- 
fore, the P and E terms in Eqs. (QJ and ((2|) can be ne- 
glected, and only the term with the large Qm + <Zl£ — 4-70 
is kept, leading to a small f=s 5% contribution [12j |. 

At order Q 3 , 2b currents enter in chiral EFT [5j. These 



include vector spatial, axial temporal, and axial spatial 
parts [15j . The first two are odd under parity, and there- 
fore can be neglected. Consequently, for the cases studied 
here, the dominant weak 2b currents only have an axial 



spatial component 
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where t~ = (j\ x t 2 )~ and the same for <r x , k = ^(p 2 — 
p 2 -p' 1 + pi) and q = ^ (pi+p' x -p 2 — p 2 ). Equation fl3J) 
includes contributions from the one-pion-exchange 03,04 
parts and from the short-range couplings d\ , d 2 , where 
due to the Pauli principle only the combination di+2d 2 = 
c D /(g A A x ) enters [with A x = 700 MeV]. 

We study the impact of chiral 2b currents in nuclei 
at the normal-ordered lb level by summing the second 
nucleon over occupied states in a spin and isospin sym- 
metric reference state or core: J° 2b = 2~2j(^ ~ Pij)^ij> 
where Py is the exchange operator. The normal-ordered 
lb level is expected to be a very good approximation in 
medium-mass and heavy nuclei, because of phase space 
arguments for normal Fermi systems at low energies [1 61 ) . 
This approximation has also been explored for chiral 2b 
currents in nuclear matter [13] , but limited to long wave- 
lengths and without connecting 2b currents and nuclear 
forces. Taking a Fermi gas approximation for the core 
and neglecting tensor-like terms (k • <rk — ^k 2 er)T~, we 
obtain the normal-ordered lb current: 
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where p = 2/c f /(37t 2 ) is the density of the reference state, 
fcp the corresponding Fermi momentum, and I(p, P) is 
due to the summation in the exchange term, 
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The contributions from 2b currents depend on the den- 
sity of the reference state. We consider a typical range 
for nuclei p = 0.10. ..0.12 fm -3 , which corresponds to av- 
erage nucleon momenta Q ~ 150 — 200 MeV, so that 
2b currents are expected to be more important com- 
pared to very light nuclei. For these densities, we have 
I(p,P) = 0.64. ..0.66, using the Fermi-gas mean-value 
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TABLE I: Quenching q of the p = GT operator pre- 
dicted by the long-range parts of 2b currents (cd = 0) for 
given C3,C4 couplings. Dots correspond to a density range 
p = 0.10. ..0.12 fm -3 , and all cz,c& values are in GeV -1 . In 
addition, we give the short-range coupling cd obtained by re- 
quiring that 2b currents lead to a quenching of q = 0.74 or 
q = 0.96 over this density range. For given q, the * values 
correspond to the weakest/strongest p dependence in Fig. [3] 
(smallest/largest | C3 |p) , which yield the thin bars in Fig. [2] 

P 2 = 6fc|/5, and I{p,P) = 0.58. ..0.60 for P = 0, so 
the total-momentum dependence is very weak. 

The effective lb current 3°% h only contributes to the 
GT operator and can be included as a correction to the 
gA (p 2 ) part of the lb current, Eq. ((2]). This result is 
general for a spin/isospin symmetric reference state. Be- 
yond the Fermi-gas evaluation, only the momentum de- 
pendence will be replaced by a state/orbital dependence. 

For cd = 0, the leading 2b currents are fully deter- 
mined by the couplings 03,04 from nuclear forces. We 
take c 3 ,c 4 from the N 3 LO NN potentials of Ref. [HI 
(EM) or Ref. (EGM), as well as from a NN par- 
tial wave analysis (PWA) extraction [2(J. Because the 
C3,C4 are large in chiral EFT without explicit Deltas, 
we also explore typical changes expected at higher or- 
der 8c 3 = — <5c 4 w 1 GeV -1 [4j. The different C3,c 4 are 
listed in Table HI As shown, the corresponding long-range 
parts of chiral 2b currents predict a quenching of the GT 
operator at p = ranging from q = 0.85. ..0.66 for den- 
sities p = 0.10...0.12fm -3 . The largest quenching is ob- 
tained from the largest 2C4-C3 coupling, which is closest 
to the single-Delta excitation with 2c 4 -C3 ~ 15 GeV -1 . 
The density dependence of the 2b-current contributions 
for cd — is also shown as inset in Fig. |3J 

This demonstrates that chiral 2b currents naturally 
contribute to the quenching of GT transitions. A re- 
duction of gA in the currents is also expected consider- 
ing chiral 3N forces as density-dependent two-body in- 
teractions 21]. We can constrain the short-range cou- 
pling cd by requiring to reproduce an empirical quench- 
ing q — 0.74 needed in many-body calculations [2]. This 
leads to values \cd\ < 1 in Table HI which is compatible 
with independent determinations from 3N forces, espe- 
cially in the best studied EM case @, [22[ (the comparison 
should be made for large 3N cutoffs); e.g., cd's obtained 
from the 3 H half-life fit @ lead to q = 0.75. ..0.67. Be- 
cause part of the quenching may be due to truncations in 
the many-body basis, we can ask what values of cd are 
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FIG. 3: (Co[or onhne) Top panel: gA plus 2b-current con- 
tributions for p — GT transitions normalized to gA as a 
function of density p. The boundaries of the box are given by 
q = 0.74/0.96 and p = 0.10. ..0.12 fm" 3 . The different curves 
correspond to the couplings in Table [T] with shaded regions 
for the density range. Intermediate q values would lie between 
these regions. The inset shows the quenching predicted by the 
long-range parts of 2b currents only (cd = 0). Bottom panel: 
Same as top, but as a function of momentum transfer p for 
empirical/smaller quenching q = 0.74/0.96. 

needed to reproduce a smaller quenching q = 0.96 (based 
on an extraction q 2 = 0.92±0.11 from GT strength func- 
tions to high energies 23j). As shown in Table HI the 
resulting cd are more negative, which seems compatible 
with 3N force fits only for the EGM and PWA cases [22] . 
Because of these uncertainties, we consider a conservative 
range from empirical to small quenching. The resulting 
quenching of gA is shown in Fig. [3] versus density. Chiral 
2b currents will universally affect p ss GT lifetimes and 
strength functions by q 2 and 2v(3(3 lifetimes by q 4 . 

The momentum-transfer dependence is predicted in 
chiral EFT and shown in Fig.[3]for the 2b-current contri- 
butions. This p dependence is determined by the density 
and the C3 coupling in Eq. (Q}, and is stronger for larger 
\cs\p values. With increasing p ~ m^, the impact of 2b 
currents, and thus the quenching, is reduced and the GT 
lb operator could even be somewhat enhanced in the 
smaller quenching case. A reduction is consistent with 
muon capture studies, which probe p y «V ~ 100 MeV, 



where GT quenching is not needed 24 1. 

Finally, we apply chiral EFT currents to the OvBj: 
cay NME, which is given by the nuclear currents [l( 
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Here, a size scale R = 1.2A 1 / 3 fm is introduced so that 
M 0v PP is dimensionless, and the closure approximation 
was used with an average intermediate-state energy E\. 
Constructing the 0^/3/3 decay operator in chiral EFT, 
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the nuclear currents are given by the lb currents J 4 ° lb and 
Jijib of Eqs. ([IJ and (J2J, and we include 2b currents at the 
normal-ordered lb level of Eq. @. With these, the 
transition operator can be decomposed into GT, Fermi 
(F), and tensor (T) terms that receive contributions from 
the V, A, P, M components in the currents [Jj| : 



Jf J 2fl = h GT (p 2 ) (7\ • o~2 h* v (p 2 ) - h T S 12 (p) , (7) 
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At leading order Q°, only the p 



vector- vector terms h^(0) and hy V \ 



- axial-axial and 
L AA\ yj I "wW contribute. At 
order Q 2 , this includes the p 2 terms in gv,A{p 2 ) and the 
other terms in the lb currents. The Q 2 terms are also 
important in phenomenological currents, as they lead to 
short-range correlation effects being small [25J, which is 
consistent with chiral EFT. At order Q 3 , 2b-current con- 
tributions enter in and h G J,. The h T parts are small 
(less than 1%, except in 48 Ca) [JJ, 14 1, consistent with 
neglecting (k • er k — ifc 2 <x) terms from the 2b currents. 

We then calculate M 0vf313 based on chiral EFT cur- 
rents for 48 Ca, 76 Ge, 82 Se, 124 Sn, 130 Te and 136 Xe. These 
are challenging calculations and require approximations. 
Here we perform state-of-the-art SM calculations using 
the coupled code NATHAN j2fj, ideal for 0+ states. The 
details of the SM spaces and interactions are given in 
Ref. [12]. Because 2b currents are independent of the 
many-body spaces, they should be included in all calcu- 
lations, and we focus on relative changes. 

Our results are shown in Fig. [2] for successive orders in 
chiral EFT. At order Q 2 , the NMEs are similar to calcu- 
lations based on phenomenological lb currents only, with 
differences mainly due to the Q 2 terms compared to phe- 
nomenological dipole form factors, as well as different A a 
values used. At order Q 3 , the contributions from chiral 
2b currents range from a 10% increase to a 35% reduc- 
tion of M 0v W. This range is obtained for the couplings 
of Table U and overlaps with the prediction based on the 
long-range parts of 2b currents. The changes at order Q 3 
agree with the p ~ 100 MeV estimate from Fig. [3] The 
lower M 0u PP values correspond to the quenching of gA 
in SM calculations of p ~ single-/3 and 2v/3f3 decays, 
but the p dependence weakens the quenching compared 
to the naive GT q 2 w 0.74 2 « 0.55 or 45% reduction. 
Finally, we note that in Delta-full EFT, the long-range 
single-Delta part of 2b currents would contribute at order 
Q 2 , which would shift the Q 2 results more into the Q 3 
range, improving the order-by-order expansion in Fig. [2] 
This presents the first study of chiral 2b currents for 
GT transitions and 0vf3f3 decay in medium-mass nuclei. 
Compared to light nuclei, their contributions are ampli- 
fied because of the larger nucleon momenta. For a spin 
and isospin symmetric reference state, the leading ax- 
ial 2b currents contribute only to the GT operator. A 
quenching of p ~ GT transitions is predicted by the 
long-range parts of 2b currents, and we have estimated 



the theoretical uncertainty from the short-range part as- 
suming an empirical to small quenching range. 

Chiral EFT predicts the momentum-transfer depen- 
dence of 2b currents that is probed in 0vj3f3 decay. This 
first calculation of the Ov/3/3 decay operator based on 
chiral EFT currents shows that 2b contributions to the 
NMEs range from —35% to 10% and should be included 
in all calculations. Future directions are: going beyond 
the lb approximation, investigating other electroweak 
transitions, developments of consistent interactions and 
operators (see also [21]), and astrophysics explorations. 
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